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Connectivity of Random 1 -Dimensional Networks 

V. Kurlin and L. Mihaylova 



Abstract — An important problem in wireless sensor networks 
is to find the minimal number of randomly deployed sensors 
making a network connected with a given probability. In practice 
sensors are often deployed one by one along a trajectory of a ve- 
hicle, so it is natural to assume that arbitrary probability density 
functions of distances between successive sensors in a segment 
are given. The paper computes the probability of connectivity 
and coverage of 1-dimensional networks and gives estimates for 
a minimal number of sensors for important distributions. 

Index Terms — Sensor networks, connectivity, probability, arbi- 
trary distribution, convolution, Laplace transform. 



I. Introduction 

Recently the problems of connectivity and coverage in 
wireless sensor networks have been extensively investigated 
[1]. One-dimensional networks are theoretically simple, but 
can be used in many practical problems such as monitoring 
of roads, rivers, coasts and boundaries of restricted areas. 
Networks distributed along straight paths can provide nearly 
the same information about moving objects as 2-dimensional 
networks, but require less sensors and have a lower cost. 

We derive the probability of connectivity of a 1-dimensional 
network containing finitely many sensors deployed according 
to arbitrary densities in contrast to [2]. We found an exact for- 
mula in the general case and explicit estimates for a minimal 
number of sensors for classical distributions. The main novelty 
is the universal approach to computing the probability of 
connectivity, which leads to closed expressions for piecewise 
constant densities approximating an arbitrary density. The 
feasibility of the proposed approach is demonstrated over 
different scenarios. We deal with densities of distances be- 
tween successive sensors, not with the distributions of sensors 
themselves, because sensors of 1-dimensional networks are 
often deployed one by one along a trajectory of a vehicle. 

Suppose that a sink node at the origin xq — collects 
some information from other sensors. Let L be the length 
of a segment, where n sensors having a transmission radius 
R are deployed. The sensor positions are supposed to be in 
increasing order, i.e. = xq < xi < ■ ■ ■ < Xn < L. Let fi{s) 
be the probability density function of the i-th distance yi — 
Xi — Xi^i. The probability that yi e [0, ^] can be computed 

as P{Q < yi < I) — J fi{s)ds. The resulting network is 



connected if the distance yi between any successive sensors, 
including the sink node, is not greater than R. 

We assume that the distances are independently distributed. 
The densities fi depend on the practical way to deploy sensors. 
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We consider the transmission radius R as an input parameter, 
because the range of available radii is often restrictive, while 
the number of sensors can be easily controlled in practice. 

The Connectivity Problem. Find the minimal number of 
randomly deployed sensors in [0, L] such that the resulting 
network is connected with a given probability. 

The Coverage Problem. Find the minimal number of ran- 
domly deployed sensors such that the network is connected 
and covers the segment [0, L] with a given probability. 

The example below shows that connectivity of networks in 
dimensions 1 and 2 are closely related. Distributing sensors 
from a vehicle along a path in a forest can result in a 
network located in a narrow road of some width W, see 
Fig. L Assuming that W < R and denoting the 2-dimensional 
positions of the sensors by (a;i, Zi), . . . , z„), where the 
n sensors are ordered by their x-coordinates, the coordinate 
Zi £ [~W/2,W/2] can be represented as a deviation of the 
i-th sensor from the central horizontal segment [0, L]. 
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Fig. 1. A 2-dimensional network distributed in a narrow road 

If the 2-dimensional network is connected, i.e. each dis- 
tance is not greater than the transmission radius R, then the 
Pythagoras theorem implies that Xi — Xi-i < \J R? — W'^ 
since \zi — Zi-i\ < W. If the 1-dimensional network of the 
sensors xi, . . . ,Xn projected to the horizontal segment [0, L] 
is connected for the new transmission radius \/R? ~ W^, then 
the original 2-dimensional network is also connected. 

Similarly, if the 1-dimensional network of projections covers 
[0, L], then the original 2-dimensional network covers the 
whole road [0, L] x [~W/2, W/2]. So if the width W of the 
road can be assumed to be less than the original transmission 
radius R, then the connectivity and coverage problems are 
reduced to the simpler problems for 1-dimensional networks. 

The paper is organised as follows. Related results on con- 
nectivity are reviewed in section II. In section III we state the 
main theorems computing the probabilities of connectivity and 
coverage. Sections IV, V, VI are devoted to explicit estimates 
of the minimal number of sensors for a uniform distribution, 
constant density with 2 parameters, truncated exponential and 
normal distribution. Appendices A-D contain proofs of the 
main theorems and corollaries including a method for com- 
puting the probability of connectivity for piecewise constant 
densities approximating any density in practice. 
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II. Related Results on Connectivity 

Many results on connectivity are asymptotic in the number 
of sensors, see [3], [4] for 2-dimensional networks. The net- 
work of n sensors in the unit disk is connected with probability 
1 if and only if the transmission radius R is proportional to 
y/ (In n)/n as n ^ oo [5], where In means the logarithm to 
the base e. These asymptotic results cannot be applied to real 
networks, because the rate of convergence is not clear. 

The standard assumption for finite networks is the uniform 
distribution of sensors. The authors of [6] suppose that sensors 
are exponentially distributed in a segment. Papers [5] and [7] 
consider sensors having the Poisson and exponential distribu- 
tion in square [0, 1]^, respectively, see also [8], [9]. 

An explicit analytical result on connectivity of finite net- 
works was obtained in [2], where n sensors are uniformly dis- 
tributed in [0, L]. In this case the probability of connectivity 
of the network was computed assuming ("7^) = for i > n 

KL/R . _ , 

i=o V « / 

The upper bound i < L/R implies that 1 — iR/L > 0, but 
the alternating inequality P'^ > is still highly non-trivial and 
can hardly be proved by combinatorial methods. This approach 
was generalised to the exponential distribution [10]. 
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Fig. 2. The probability of connectivity for 2 uniformly distributed sensors 

By the formula above for n = 2 sensors having a trans- 
mission radius R, the probability of connectivity is = 
1 - (1 - R/Lf = 2{R/L) - {R/Lf. This is illustrated in 
Fig. 2, where a network of 2 sensors at xi,x2 is represented 
by a point in the triangle {0 < xi < X2 < i}. Then P2 is the 
area of the domain {0 < a;2 — xi < i?} fl [0, L]^ of coimected 
networks divided by the area of the triangle. 

III. New Theoretic Results 

Recall that one deploys n sensors having a transmission 
radius R in [0, L] in such a way that the i-th distance Xi — Xi-\ 
between successive sensors has a probability density function 
fi{s) for i ~ 1, . . . , n. Assume that the densities /i, . . . , /„ 

L 

are integrable and / fi{s)ds = 1, z = 1, . . . , n. Hence the z-th 




distance can take values from to L. So the n-th sensor may 
not be within [0, L] and its position is bounded only by nL. 
A network is proper if all sensors are deployed in [0, L\. In 
practice all networks are proper, because sensors are deployed 
along a fixed segment. A proper network is connected if the 
distance between any successive sensors, including the sink 
node at 0, is not greater than R, see Fig. 2. 

We will compute the conditional probability that a proper 
network is connected, i.e. the probability that the network 
is connected assuming that it is proper. So the answer will 
be a fraction, the probability that the network is proper and 
connected over the probability that the network is proper. The 
numerator and denominator will be evaluations of the function 
Vn{r, I) defined recursively for n > as follows: 

wo(r, = 1 if r, Z > 0; 

Vn{r,l)=Q ifr<OorZ<0; 
Vn{r,l) = l ifr>Z>0, n>0; 

r 

Vn{r, I) = J fn{s)vn-i{r, I - s)ds if r < I, n > 0. 


The Probability Proposition. For < r < Z in the above 
notations, Vn{r,l) is the probability that an array of random 
distances (yi, . . . , y„) with densities fi,. . . ,fn, respectively, 

n 

satisfies J2 Vi — ^ < Ui < r for i = 1, ... ,n. 

i=l 

The variables r, I play the roles of the upper bounds for the 
distance between successive sensors and the sum of distances, 
respectively. Clearly 77„(L, L) is the probability that a network 
is proper, i.e. all sensors are in [0,L], and Vn{R,L) is the 
probability that a network is proper and connected. 

The Connectivity Theorem. Let n sensors xi , . . . , a;„ having 
a transmission radius R be deployed in [0, L] so that a sink 
node is fixed at xq = and the distances yi = Xi — Xi-i, i = 
1, . . . , n, have given probabiHty density functions /i, . . . , /„. 
Then the probability of connectivity of the resulting network is 

Pn = ^" ^ which is independent of the order of sensors, 

Vn(L,L) 

the function Vn{r, I) was recursively defined above. 

Given a probability p, the answer to the Connectivity Prob- 
lem from section I is the minimal number n such that P„ > p. 
A network of a sink node at and 1 sensor with at i/i e [0, L] 
is coimected with probability Pi = P(0 < yi < R) = 

R L 

= vi{R,L) = Jfi{l)dl, since vi{L,L) = j f^{l)dl = 1. 



The Coverage Theorem. Under the conditions of the Connec- 
tivity Theorem, the probability that the network is connected 

, wnn- v„{R, L) - v„{R, L - R) 
and covers the segment 0, L\ is , ^ ^, . 

v„(L,L) 

The Connectivity Theorem leads to closed expressions for 
probability of connectivity and explicit estimates on a minimal 
number of sensors making a network connected with a given 
probability for classical densities in sections IV- VI. The Con- 
nectivity and Coverage Theorems are proved in Appendix A 
by generalising the analytical method from [2]. 
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IV. The Uniform Distribution 

In this section we consider the simplest constant density 
f{l) = 1/L on [0, L], i.e. the distances between successive 
sensors are uniformly distributed in [0, L]. The formula for P," 
in the Uniform Corollary below can be compared with the for- 
mula for from section II obtained in [2] for networks whose 
sensors (not distances) are uniformly distributed in [0,L]. In 
the latter case there is no sink node at 0, see differences in 
Fig. 2-3. In Fig. 2 the network of 2 sensors is represented by 
thek positions {xi,X2), while in Fig. 3 the same network is 
encoded by the distances (2/1,2/2) — {xi — 0,X2 — xi). 

The Uniform Corollary. Under the conditions of the Connec- 
tivity Theorem, if the distances between successive sensors are 
uniformly distributed in [0,L], then the probability of connec- 

tivity is P,^ = '^e''(-1)'C;)(1 - tR/ir. Set = ^ - 1. 
The network is connected with a given probability p > 2/3 if 




Fig. 3. The probability of connectivity for a sink node at and 2 sensors 
with uniformly distributed distances yi = a;i — and 3/2 = X2~xi in [0, L] 

For n = 1 the Uniform Corollary gives P" = R/L, namely 
a network of a sink node at and another n — 1 sensor at a 
distance j/i = — is connected if and only if j/i < R, i.e. 
with probability = R/L. For n = 2 one gets: 

r 2(P/L)2 if P < L/2, 

2 - \ 4(P/P) - 2(P/L)2 - 1 if P > L/2. 

If R < L/2, then the probability is the area of the square 
{0 < yi < P, < y2 < R} divided by the area of the 
triangle {0 < yi, < y2, < j/i + ?/2 < L}, see Fig. 3. The 
lower bound in the Uniform Corollary is positive if P > 2P, 
because the 2nd term under the square root is non-negative 
for p e (0, 1) and the square root is not less than 3Q — 1. 

The computational complexity of P„ is linear in the number 
n of sensors. By the computational complexity we mean 
the number of standard operations like multiplications and 
evaluating simple functions like ln(a;). A linear algorithm 



computing P„ above initialises the array consisting of n + 1 
elements L — iR, i ^ 0, . . . ,n, then finds ln(i — iP), 
n ln(P — iR) and exp(n ln(P — iR)) = (i — iP)". The array 
of binomial coefficients (") has n + 1 elements and can be 
computed in advance. So the total complexity of computing 
the probability in the Uniform Corollary is 0{n). 

Consider the segment of length L = 1km and n < 200 
sensors having transmission radius P = 50m. Suppose that a 
sink node is fixed at and the distances between successive 
sensors have the same uniform distribution on [0, P]. The 
graph in Fig. 4 shows the probability P^ of connectivity 
computed in the Uniform Corollary. The number n of sensors 
varies from 1 to 200 on the horizontal axis. 
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Fig. 4. The probability of connectivity for the uniform distribution 

The graph in Fig. 4 implies that after a certain value of 
n the probability P^ of connectivity increases with respect 
to the number of sensors. To solve the Connectivity Problem 
from section I for a given probability p, we compute P^ for 
all values from 1 to a minimum n such that P^ > p. Another 
method uses the estimate from the Uniform Corollary, which 
may not be optimal, but requires much less computations. The 
exact minimal numbers and their estimates are in Table 1, 
where the network in [0, L] with L — 1km is connected with 
probability p = 0.95. For example, the minimal number of 
sensors for P = 50m is 157, while the estimate is 905. 



Table 1. Simulations for the uniform distribution 



Transmission Radius, m. 


200 


100 


50 


25 


10 


Min Number of Sensors 


29 


69 


157 


349 


982 


Estimate of Min Number 


83 


283 


905 


2610 


10640 



Table 1 implies that the uniform distribution is very ide- 
alised and can not be useful in practice. If one deploys sensors 
of transmission radius P = 50m non-randomly at regular 
intervals 49m, than 2 1 sensors are enough to make the network 
connected and cover [0, P]. The estimate from the Uniform 
Corollary is too rough, because of the term 24(5'^/(l — p) 
under the square root, which can be very large. The uniform 
distribution is extended to a more practical case in section V. 
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V. A CONSTANT DENSITY WITH 2 PARAMETERS 

In this section we consider the constant density f{l) = 
l/(&— a) over any segment [a,b] C [0, L], which generalises 
the uniform distribution from section IV. Practically the distri- 
bution means that each sensor is thrown at a distance uniformly 
varying between a and b from the previously deployed sensor 



f{l) 



a 



1/(6 -a) 




R 


b 



Fig. 5. The constant density over the segment [a, b] 

The left endpoint a should be less than the transmission 
radius R, otherwise no sensor communicates with its neigh- 
bours. The mathematical expectation of the distance between 
successive sensors is (a + h)/2, while (6 — a)^/12 is the 
variance of the distance. For example, for a network of a sink 
node at and 1 sensor at j/i, the probability of connectivity 
is P(0 <yi<R)^{R- a)/{b - a), see Fig. 5. 

The Constant Corollary. If in the Connectivity Theorem the 
distances between successive sensors have the density f{l) — 
l/{b — a) on [a, b], then the probability of connectivity is 



E(-l)'^a)(i-a(n-fc)-i?fcr 
fc=o 



E("i) 

k=0 



I'd) {L - a{n - k) - bky 



The network is connected with a given probability p if 



n > max < — 
2 



1 + 5p 



1 



L-b 



and 



< R<b. 



The sums include all expressions taken to the power n if 
they are positive. The complexity to compute is 0{n). 
Each of the terms in both sums requires 0(1) operations 
similarly to the Uniform Corollary. For n = 1 one gets 

Pf . = ^ as expected above. If 

[L — a) — [L — bj b — a 
the given probability p is too close to 1 then the estimate 

from the Constant Corollary depends on p, e.g. n > 776 

for p = 0.9999, but in all reasonable cases the maximum is 

achieved at the second expression l + {L — b)/a independent of 

a + b 

p. The restnctions — - — < R < b seem to be natural saying 
that the distance between successive sensors is likely to be 
less than R since [0, R] covers more than a half of [a, b]. 

Each distance Xi — Xi-i between successive sensors belongs 
to [a,b]. Such a network lies within [0,L] only if an < L, 
hence the number of sensors should satisfy n < L/a. In 
the boundary case an = L all sensors should be located at 
the exact positions Xi — ia/n, i ~ 1, . . . ,n, which clearly 



happens with probability 0, so the numerator vanishes for 
L = an in the Constant Corollary. If n > L/a then n sensors 
can not be within [0, L] according to the density 1/(6 — a). 




Fig. 6. The probability of connectivity, R = 50m, a = 0.2R, b 

Table 2. The case of the constant density over 



Transmission Radius, m. 


200 


150 


100 


50 


25 


Min Number of Sensors 


14 


19 


30 


63 


132 


Estimate of Min Number 


18 


27 


43 


93 


193 


Max Number of Sensors 


25 


34 


50 


100 


200 



: i.m 



0.2P, 1.6P1, 



Figs. 6-8 show the probability of connectivity for different 
segments [a, b] depending on the radius R = 50m. The graph 
in Fig. 6 is the probability of connectivity for 1 < n < 100, 
i 1 km, P = 50m, a = 0.2P, b = 1.6P. If the required 
probability of connectivity is p = 0.95 and the transmission 
radius is 50m, then the minimal number of sensors is 63. 

The maximal possible number of sensors is L/a = 100, i.e. 
P^QQ = since the sensors should be fixed at exact positions 
in [0,L], which explains the drop to in Fig. 6. The minimal 
number of sensors decreases when the length 6 — a decreases. 
The maximum number of sensors in Table 2 is L/a, which 
gives probability in this extreme case. All numbers slightly 
less than the maximum give a probability close to 1. More 
exactly we may subtract b/a ~ 1 = 7, see Table 2, which 
follows from the second restriction n > 1 + {L — b) / a. 

Table 3. The case of the constant density over [0.4P, 1.4P]. 



Transmission Radius, m. 


200 


150 


100 


50 


25 


Min Number of Sensors 


10 


13 


20 


41 


83 


Estimate of Min Number 


10 


14 


22 


47 


97 


Max Number of sensors 


13 


17 


25 


50 


100 



Table 4. The case of the constant density over [0.6P, 1.2P] 



Transmission Radius, m. 


200 


150 


100 


50 


25 


Min Number of Sensors 


8 


10 


15 


31 


61 


Estimate of Min Number 


8 


11 


15 


32 


65 


Estimate of Max Number 


10 


13 


17 


34 


67 



Tables 2-4 imply that the required number of sensors 
making a network connected decreases if the ratio {b — a)/R 
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Fig. 7. The probability of connectivity, R = 50m, a = 0.4R, b = 1.4R 
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1 < n <33 



Fig. 8. The probability of connectivity, R = 50m, a = 0.6R, p = 1.2R 

decreases. For b — a < R the estimate from the Constant 
Corollary is very close to the exact minimal number of sensors 
when sensors are deployed non-randomly at a distance slightly 
less than R. So the found estimate for the minimal number of 
sensors requires few computations and can be useful. 

VI. Exponential and Normal Distributions 

Here we state partial results for 2 other classical distribu- 
tions. For the exponential density over [0, L], we compute the 
exact probability of connectivity, but the simplest estimate for 
the minimal number of sensors is the same as for the uniform 
distribution. For the normal density, it is hard to compute the 
probability of connectivity explicitly, but a reliable estimate 
for a maximal number of sensors will be derived. 

Consider the exponential distribution /(s) = Ae^'*'*, A > 0. 
It is used for modelling the wait-time until the next event in 
a queue. Since sensors are deployed in [0,L], we consider 
the truncated density f{s) = ce^'*''' on [0, L] and f{s) — 

A 



The Exponential Corollary. If in the Connectivity Theorem 
the distances between successive sensors have the exponential 
density /(s) — ce^^" in [0, L], then the probability of 



connectivity is 



Vn{R, L 
Vn{L, L] 



, where Vn{r, I) 



i<l/r 

E 

i=0 



(-1)' 



-iXr 



A" 



E 

j=0 



The estimate for a minimal number of sensors from the 
Uniform Corollary holds in this case, which can be proved an- 
alytically, but easily follows from the fact that the exponential 
density monotonically decreases on [0, L], hence the distance 
between successive sensors will be smaller on average than 
for the constant density over [0,L], i.e. the network is more 
likely to be connected. The computational complexity of in 
Corollary 2 is O(n^), because each expression in the brackets 
requires 0{n) operations as in Corollary 1 assuming that ln(.T) 
and exp(a;) can be computed in 0(1) operations. 

The Exponential Corollary implies that u„(L,L) = 1 — 

n-l 

g-AL [XLy / j\ Indeed the term corresponding to i = 1 

i=o 

1 



vanishes if r = 1. 



The sum ^ 



converges rapidly 



to e^'^'"*'"^ as n ^ oo. Hence the expression in the brackets 
from the Exponential Corollary is very close to even for 
small n. Then is a ratio of tiny positive values of order 
10~^° or less. The computation of P^ very fast accumulates 
a big arithmetic error even for small n. The exponential 
decreasing of ce~^^ means that the sensors are distributed very 
close to each other and cover [0, L] with little probability. So 
the exponential distribution seems to be rather unpractical for 
modelling distances between successive sensors. 

Finally we consider the remaining classical distribution, 
the truncated normal density over [0,-L], i.e. /(s) = 
/2o- ^ where the constant c guarantees that 



(TV27r 

/o f{s)ds — 1. The normal density has exponentially decreas- 
ing tails, so distances between successive sensors are likely 
to be close to /i. Hence the mean /i should be less than the 
transmission radius R and the number of sensors n can not 
be greater than L/fi, otherwise last sensors are likely to be 
outside [0, L]. That is why the Normal Corollary below gives 
an upper bound for the number of sensors making a network 
connected, not a lower bound as in previous corollaries. 

The Normal Corollary. If in the Connectivity Theorem the 
distances between successive sensors have the truncated nor- 
mal distribution on [0, L] with a mean /i and standard deviation 
a then the network is connected with a given probability p for 



n < min ■ 



p{l~p) (V4^L + c72$-2(p)-a$-i(p))^ 
e ' 4/i2 



otherwise. The condition J f{s)ds ~ 1 gives c 



1 



6 



The standard normal distribution $(a;) is not elementary, 
but its values have been tabulated. The table below shows 
estimates for the maximal number of sensors normally dis- 
tributed in [0,L] with L = I km, (i = 0.6R, a = O.IR 
in such a way that the resulting network is connected with 
probability p = 0.9975. Then w 2.8, e w 0.000063 

and the first upper bound in the Normal Corollary gives 
n < p{l — p)/e « 40, which is the overall upper bound 
for R = 25m. For radii R > 50m the second upper bound is 
smaller that the first one and is close to L/ ji, the exact number 
of sensors when all distances are not random and equal to /x, 
because a^~^{p)/R ~ 0.28 is rather small. 

Table 5. The case of the normal density, fi = 0.6R, a = O.IR. 



Transnoission Radius, m. 


200 


150 


100 


50 


25 


Estimate of Max Number 


7 


11 


16 


33 


40 



The estimates from Table 5 are close to optimal, e.g. for 
the radius R = 150m the non-random distribution of sensors 
at distance 149m apart requires 6 sensors not including the 
sink node at 0, while the estimate above gives 11. The ratio 
6/11 is close to the mean ij,/R = 0.6 since distances between 
successive sensors should be around the average /j. = 0.6R. 

VII. Conclusions 

We would like to emphasise that the main result of the 
paper is a new method of analytical computing the probability 
of connectivity of random 1 -dimensional networks leading to 
explicit formulae for piecewise constant densities approximat- 
ing an arbitrary density. The found estimates for a minimal 
number of sensors making a network connected suggest that 
a constant and normal densities over a segment can be more 
economic than other other classical distributions. 

Open issues for the future research are the following: 

i) computing analytically the exact probability of connectivity 

in the case when the distances between successive sensors have 
a truncated normal distribution over [0,L]; 

ii) finding an optimal distribution of distances between suc- 
cessive sensors in [0,1/] for a given number of sensors to 
maximise the probabilities of connectivity and coverage; 

Hi) extending the suggested approach of sensor distributions 
to non-straight trajectories filling a 2-dimensional area. 

Appendix A 
Proofs of the Main Theorems 

First we recall the notion of the convolution and Laplace 
transform used in the proof of the Connectivity Theorem, 
Coverage Theorem and Corollaries from sections IV- VI. The 

+ 00 

convolution of functions /, g is / * g{s) = / f{l)!){s — l)dl- 

— DC 

The convolution is commutative, associative, distributive and 
respects constant factors, i.e. {cf)*g = c{f*g), f*g = g*f, 
{f*g)*h = f*{g*h), f*{g+h) = f*g+f*h. The convolution 
plays a very important role in probability theory, because the 
probability density of the sum of 2 random variables is the 
convolution of the densities of the variables. 



Given a function f{l) and r > 0, introduce the truncated 
function = f{l) for I € [0, r] and = otherwise. 

Let u{l) be the unit step function equal to 1 for /> and 
equal to for Z < 0. Then the truncated function f^'~\l) is 
f{l){u{l) — u{l — r j). Below we use the partial convolution 
/(r, I) * g{r, I) considered only for the argument I, while r 
remains constant. The following lemma rephrases the recursive 
definition of Vn(r,l) in terms of convolutions. 

Lemma 1. Given densities /i,...,/„, the function Vn{r,l) 
from Section HI is fn^ * ... * /j'' * u{l), r < I, n > 0. 

Proof of Lemma 1 is by induction on n. The base n = 1 is 

trivial: * u{l) = J /i(s)u(s — l)ds = J fi{s)ds = vi{r, I) 



since s <r < I. The inductive step follows from the recursive 
definition of w„ in section III: w„(r, I) = fn^ * Vn-i{r, I). □ 

The Laplace transform of a function f{l) is the function 

LT{/(/)}(s) = / e-'^f{l)d.l. The Laplace transform is a 



linear operator converting the convolution into the product, 
i.e. LT{a/ + bg} = aLT{/} + bhi:{g}, LT{/ * g} = 
LT{/}LT{g'}. The inverse Laplace transform LT ^ is also 
a linear operator. The following well-known properties of the 
Laplace transform can be easily checked by integration. 

Lemma 2. For any a, f3 and integer m > one has 

(a) LT{l-u{l)} = 

T 7? ^ 

(b) LT{e-"'/™u(/)} - '■ 

[s + a)™+^ 

^lp-/3s 

(c) LT{(/ - pru{l - /?)} = and 

(d) LT{e-('-/5)(Z - (3ru{l - (3)} = j^^^^- □ 

Lemma 2 allows one to compute the inverse Laplace 
transform, e.g. Lemma 2(a) implies that LT~^{l/,s} = u{l). 
Lemma 3 provides a powerful method for computing the 
function Vn{r, I) used in the Connectivity Theorem. 

Lemma 3. Given probability densities /i, . . . , /„ on [0,L], 
set g{s) = LT{/Jr'(0} ' • • • ' LT{/fl(0}/s. Then the function 
Vn{r,l) is the inverse Laplace transform LT^^{g(.s)}(Z). 

Proof of Lemma 3. One has t'„(r, I) = fn^ *•••*/-[''* u{l) by 
Lemma 1. Set g{s) = LT{w„(r,Z)} and g,{s) = LT{/fl(0}, 
i = l,...,n. The Laplace transform is considered with 
respect to I, the variable r is a fixed parameter. The Laplace 
transform converts the convolution into the product, hence 
g{s) = gi{s) . . . gn{s) / s as expected since LT{u{l)} = 1/s 
By Lemma 2(a). The order in the product g{s) does not 
matter as the convolution is commutative. So any reordering 
of the densities gives the same result and the probability of 
connectivity does not depend on this order. □ 

Proof of the Connectivity Theorem. 

Let = xo < a;i < • • • < x„ < L be the positions 
of a sink node and n sensors. Suppose that the distances 
yi = Xi — Xi-i, i = l,...,n, are independent and have 
probability densities fi{s). Any network can be represented by 
ordered sensors {xi, . . . ,Xn) or, equivalently, by the distances 
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{Ui, . . . , Hn) between successive sensors. Then the conditional 
probabiHty of connectivity is the probabiHty that the network 

n 

is proper and connected, i.e. J2 Hi — ^ '^^'^ < Ui < R, 

i=l 

divided by the probabiHty that the network is proper, i.e. 

n 

< L and < y-i < L. Hence the required formula 

i=l 

Vn(R, L) 

Pn = — 7 r for the conditional probability of connectivity 

Vn(L,L) 

follows from the Probablity Proposition stated in section III. 
Permuting densities leads to the same probability P„ due to 
commutativity of the convolution from Lemma 1 . □ 

Proof of the Probability Proposition. 

We illustrate the proof first in the partial cases n = 1,2. For 

r 

n = 1 and L > r, P{0 < yi < r) = J fi{s)ds = vi{r, I) and 



I 

-P(0 <yi<l) = J fi{s)ds = vi{l,l) as expected. 



For n = 2, let the distance y2 belong to [s, s + A] C [0, r] 
for some small A > 0. The probability of this event E is 
P{E) = P(s < 2/2 < s + A) « /2(s)A, the area of the 
narrow rectangle below the graph of /2 over [s, s + A]. The 
random variables yi = xi — xq and y2= x-2— xi are assumed 
to be independent. Then the probability of connectivity is 

P{E)P{0 <yi<l-s)^ /2(s)A • v^ir, I - s). 

The total probability is the limit sum of the above quantities 
over the intervals [s,s + A] covering [0, i?] when A — > 0. 

r 

Hence the probability is / f2{s)vi{l — s)ds = V2{r, I). 



We will prove the general case n > 1 by induction on n. 
If the network is proper and connected then the nth distance 
yn = Xn — Xn-1 < is not greater than r and not greater 

n-l 

than I — J2 Hi- The former condition means that the last 

i=l 

sensor is close enough to the previous one. The latter condition 
guarantees that all the sensors are in the segment [0,^]. 

Split [0, r] into equal segments of a small length A. Suppose 
for a moment that /„(s) is constant on each segment [s, s+A], 
where s — jA, j = 0, . . . , [r/A] — 1. The general case will 
be obtained by taking the limit under A — ^ 0. 

The probability P{yn & [s,s + A]) is approximately 
/„(s)A, the area below the graph of fn{s) which is assumed 
to be constant over a short segment [s, s + A]. The probability 
that the n — 1 sensors form a connected network in [0, ^ — y„] 
is approximately Vn-i{r, I — s) hy the induction hypothesis. 

Since the distances are distributed independently, the joint 
probability is /„(s)A • Vn-i{r, I — s). The total probability is 
Vn{r, I), the limit sum over aU these events as A — > 0: 

[r/A] r 

51 fnU^)^ ■ Vn-i{r, I - jA) / f„{s)v„-i{r, I - s)ds. 

i=i 

The final expression above is the standard definition of the 
Riemann integral of fn{s)vn-i{r, Z — s) as a limit sum. □ 

Proof of the Coverage Theorem. 

By the Probability Proposition (i?, L) is the probability of 



the event E{L) that n sensors are deployed in [0, L] and form 
a coimected network. The network covers [0, L] if also at least 
one sensor is in [L — R, L], i.e. E(L — R) does not happen. 
Hence the probability that the network is proper, connected 
and covers [0, L] is Vn{R, L) — ?;„(i?, L — R). So the required 
conditional probability assuming that all sensors are in [0, i] 

. Vn{R,L)-Vn{R.L-R) . 

IS equal to ^. as required. □ 

Vn{L,L) 

Appendix B 
Proofs of the Main Corollaries 

The iterated convolutions respect constant factors, i.e. 

(c„/„)M * • • • * (ci/i)M *u = cn... ci/M * • • • * 1 * u. 

Hence we may consider probability densities without extra 
factors if we are interested only in the conditional probability 
Pn from the Connectivity Theorem. Indeed, the product of 
these factors will cancel dividing Vn{R, L) by Vn{L, L). 

Proof of the Uniform Corollary. 

Let f{l) be the uniform density over [0, L], i.e. /(/) = 1/L on 
[0,L], /(/) = otherwise. Let * g be the n-th iterated 
convolution, e.g. f^^*^ * 9 = f * if * g)- 

Lemma 3 gives a straightforward method to compute 
Vn{r,l) = * u, where u{l) is the unit step function, 

i.e. u{l) = 1 for Z > and u{l) = for Z < 0. 

Assume that /['"^(Z) = u{l) — u{l — r) forgetting about front 

factors. Lennma 2(c) impUes that LT{/[''1} = . By 

s 

Lemma 3 one has Vn{r, I) = LT {g{s)}, where 
(1-e 



9{s) = 



i=0 ^ ' 



By Lemma 2(b) LT"^ 



„n+l 



Replacing u{l — ir) by the upper bound i < l/r, we get 



i<l/r 



n\ {I — ir)" 



Vnir,l) = LT-\g{s)}{l)= (-1)' 



By the Connectivity Theorem the denominator of P^ = 
Vn{R,L)/vn{L,L) is Vn{L,L) = L"/n\ Hence we may 
divide each term (L — iR)"/nl in ?'„ (-??, L) by L"' /n\, which 
gives the final formula from the Uniform Corollary. 

Now we prove the estimate for a minimal number of sensors 
making the network connected with a given probability p, i.e. 
we should check that the probabiUty > p if 



n>\ 1^3(1 -Q) + y(3Q- 1)2 + 24Q2 . 

where Q = (L/R) — 1. The idea is to simplify the inequality 
Pn > P replacing P^ by smaller and simpler expressions, 
which will lead to the required lower bound for n above. 
Setting q = R/L, the probability from the Uniform Corollary 
becomes the alternating sum starting as follows: 



P„" = l-n(l-# 



(l-2# 



(l-3(?)" + . 
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The sum involves only positive terms of the form 1 — iq. 

First we check that > 1 — n(l — g)" forgetting about 

the remaning terms. It suffices to show that every odd term 

(2fc+i)(l — {2k + l)q)^ is not greater than the previous even 

one (2k) ~ 2^9)" for /c > 1. The last inequahty is equivalent 

l-2kq Y n-2k ^ , , , . , , . , 

Replace the left hand side 



to 



> 



l-{2k + l)qj - 2k + l 
by the smaller expression (1 — and the right hand side 
by the greater expression n/3 using fc > 1. 

The resulting inequality (1 — 5)"" > n/3 is weaker than the 
simplified inequality P" > 1 — — q)" > p forp € (2/3, 1), 
i.e. > n/(l-p). We check that > n/{l-p) 

holds under the required restriction on n. Since 

Q = (L/R) - 1 = (1 - q)/q then (1 - g)"" = (1 + l/Q)". 

Therefore the proof finishes by the following lemma. 

Lemma 4. The lower bound for n from the Uniform Corollary 
implies that (1 + l/Q)" > n/(l - p). 

Proof of Lemma 4. Expand the Taylor series (1 + l/Q)" = 

= 1 + n/Q + n{n - 1)/2Q2 + n{n - l){n - 2)/6Q^ + ■■■ 

Leaving the terms of degrees 1,2,3 only makes the inequality 
(1 + l/Q)" > n/(l — p) stronger, hence it suffices to prove 

n/Q + n(n - 1)/2Q^ + n(n - l)(n - 2)/6Q^ > n/(l - p). 

Multiplying boths sides by QQ^/n, we get 

6Q2 + 3Q(n - 1) + (n - l)(n - 2) > 6QV(1 - p), 

(n - 1)2 + (3Q - l)(n - 1) + 6Q\l - Q/(l - p)) > 0. 

The quadratic inequahty holds if n — 1 is not less than the 2nd 

root ^ ^1 - 3Q + y'(3Q - 1)2 + 24Q2 (^3-3^ 

which is equivalent to the required condition on n. 

Proof of the Constant Corollary. 

The truncated constant density over [a, b] without extra factors 
is = u{l - a{r)) - u(l - b{r)), where [a(r),6(r)] = 

[a, h] n [0, r] is the domain, where the probability density is 
defined and restricted to [0, r]. For instance, if < a < i? < 
b<L then a{L) = a{R) = a and b{L) = b, b{R) = R. 

By the Connectivity Theorem and Lemma 1 the probability 

is expressed in terms of v„(r, /) = (/['"1)("*) * u{l). 
Lemma 2(c) for rn = 0, /3 = a(r), P = b{r) impHes that 




LT{u(Z - a(r)) - u{l - b{r))} = 



^-a{r)a _ ^-b{r)s 



Apply Lemma 3 multiplying n factors and dividing by s: 



9{s) = 



E(-i) 

k=0 



^g— o(r)s _ g— 5(r)s^n 



-^ — a{r)s{n — k)—b{r)sk 
„n+l 



After expanding the binom, compute the inverse Laplace 
transform of each term by Lemma 2(d) for the parameters 
a = 0, f3 = a{r){n— k) + b{r)k, m = n as follows: u„(r, /) = 

= LT->b(.)} = EC-D'Q "-°'-"-„f'-^"''"'" . 

fc=0 ^ ^ 

To get the final formula for the conditional probability P^ = 
Vn{R,L)/vn{L,L) of connectivity it remains to substitute 
a{L) = a{R) = a, b{L) = b, b{R) = R and cancel n! 

Now we prove the estimate for a minimal number of sensors 
making the network connected with a given probability p. The 
condition n > l + {L — b)/a, i.e. L — a{n — 1) — 6 < 0, implies 
that the denominator of P^ from the Constant Corollary 
is equal to {L — anY^ corresponding to rn = 0. Another 
assumption {a + b)/2 < R means that b < 2R, hence 
the numerator of P^ contains only the first 2 terms, namely 
(L — an)" — n(-L — a(n"l) — i?)". The equivalent inequalities 



L-a{n-l)-RY 



an 



>P, 



L — an 



L-a{n-l)-Rj ~ 1-p 



> 



1 + 



R-a 



L-a{n-l)-Rj ~ 1-p 



> 



n 



f _ R — a\ n ^ , . , , 

are weaker than 1 H > as L — a(n— 1) < 0. 

V h- RJ - 1-p ^ ' - 

^ 

By Lenmia 4 for Q = — the last inequality holds if 

R — a 



n>i 1^3(1 -Q) + ^(3Q -1)2 + 24Q2 . 

Since {a + b)/2 < R then Q < 1 and we replace Q by 1 in 
the last expression making the condition on n only stronger: 




l + 5p 
1-p ■ 



□ 



Proof of the Exponential Corollary. 

By the Connectivity Theorem and Lemma 1 it suffices to 

compute Vn{r,l) = (/['"1)("*) * u{l), where the probability 
density function is /''"'(O = e~^''{u{l) — u{l — r)) without 
extra factors. Lemma 2(d) for a = A, /3 = r implies that 



LT{e-^^'-''>u{l-r)}{s) = 



s + A' 



LT{/H(0}(s) = 



1 _ e-''(«+^) 
s + A ■ 



By Lennma 3 one has Vn{r, I) = 1ST ^{g{s)}, where 



9is) = 



E(-i)' 



The following result will be easily proved later. 
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Lemma 5. For any A > and n > one has 



1 I ^ I 



— -T 



s{s + A)" A"s ^ A"--? (s + A)-?+i ■ 



By Lemma 2(d) for a = A, /3 = ir one has 



J /M (s)P(^ < ; - s)ds - (n - 1)^ j > 

-ex. »=1 

+?° n-1 . 

> J /(s)P(^2/z</-s)ds- y f{s)ds-in-l)e 



IR-[0,r] 



n-1 



It remains to apply Lemma 5, collect all terms in one sum and — J/i < ~ '^^ using f 2/i < ^ — s) < 1- □ 

replace u{l — ir) by i < l/r, i.e. Vn{r, I) = LT~^{g{s)} = 

i<i/'-,_,w/x / "-'^ _ ■ \j\ Proof of the Normal Corollary. 

= ^ L_LL|'^j I 1 _ e"^('-*'') ^ — A__!!jL j . □ By the Connectivity Theorem the probability of connectivity is 

i=o \ j=o •^^ / P„ = w„(i?, (L, Z/), where the denominator i) = 

* is computed using the truncated normal 

density over [0,L], while in ^;„(i?,L) = (/[«!)("*) *w(L) 
/"roof of Lemma 5 is by mduction on n. The base n = 1 j v • * j tt, n m oi 

■' the same density is truncated over the smaller range [0,R\. 

= is absolutely trivial. As usual we may forget about extra constants in front of 

s{s + X) As A(s + A) 1 / (s_„)2\ 

f{s) = — = exp — I . For a given probabiUty p 



The induction step from n — 1 to n uses the base for n = 1: (7^/2tt \ y 

we will find a condition on n such that P„ > p. We will make 
1 / 1 !lzif 1 \ 1 the inequality P„ > p simpler and stronger replacing L) 



n-2 

i V 



g(^g _(_ X)n y A"~^s ^ X"-~^~^{s + Xy j s + X Vn{R, L) by their upper and lower bounds, respectively. 

The denominator Vn{L,L) is the iterated convolution of 
1 1 1 normal densities truncated over [0, L]. This convolution of 

~ A"s ^ A"(s + A) ~ X"-~Hs + Xy+^ ' ^ positive functions becomes greater if we integrate the same 
•'=^ functions over M. Then v„{L,L) < PCEHi 2/* < 

probability that the sum of n normal variables with the mean 
In the proof of the Normal Corollary we apply the following ^ deviation a is not greater than L. The sum E-Li Vi is 
estimate for iterated convolutions of truncated probabihty ^^e normal variable with the mean n/x and deviation a^^. 
densities using tails the normal density /(s) over R. /L \ 

1 / - Th^° ^(Si-i Vi < L) = ^ ( ^ ), where the stan- 

Lemma 6. Let f{s) = ^ exp (-%?^ j be ^^"^ V / 

the normal density with I^mean n and deviation a. Then dard normal distribution is $(,t) = / er^'/'^ds. Taking 



27r_c 

" into account the lower estimate of Vn {R, L) from Lemma 6, 

Vn{r, I) = * u{l) > yi<l)- ne, where we replace the inequality P„ > p by the stronger one 



i=l 



ne 

> 



1-p 



r /7-,/V~^ T\ ji. f ^ ^ ''^l^\ 

£ = 1 - / f{s)ds and yi have the density f{s) over R. 1 - ns/P{2_^ yi<L)>por^l ^ 1 

Proof of Lemma 6 is by induction on n. The base n = 1: gp^^ ^j^^ j^^^ inequality into two simpler ones: 

+00 -\-oo 

vi{r,l)= [ /M(sM/-s)ds= / f{s)u{l-s)ds- $ ("^l—^^ > p and p > 

J J \ (Jy/n J 1-p 

— oo — oo 

The latter inequality gives n < p{l—p)/e as expected, where 

J (ti)Uti = '±'1 — "±"1 

M-[0,il] 

since u{l—s) < 1, / f{s)u{l—s)ds < / fis)ds = e. 

R-[0,r] R-[0,r] 

The induction step from n - 1 to n is similar: L-n^i> cr*"^ (p)\AI, + cr^-^ (p) - L < 0. 

"^^^ The final condition says that n is not greater than the square 

Vn{r, l)= /M(s)w„-i(r, I - s)ds = of the 2nd root + cr2$-2(p) - (7$-i(p))/2/x. □ 



i ine latter inequality gives n ^p[i.—p)/e as expectea, 

- / f{s)u{l- s)ds> / f{s)ds-e = P{yi<l)-e If / u\ fR- 

M-[0,il] 

The former inequality above becomes the quadratic one: 

-[0,r] R-[0,r] 
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Appendix C 
Networks with sensors of different types 

We derive an explicit formula and algorithm for computing 
the probability of connectivity when distances between succes- 
sive sensors have different constant densities. These general 
settings might be helpful for heterogeneous networks contain- 
ing sensors of different types, e.g. of different transmission 
radii. Assume that each distance between successive sensors 
has one of k constant densities fj(l) = Cj on [aj,bj\ C [0, L] 
and fj{l) = otherwise, j = l,...,k. The condition 

L 

J fj{l)dl = 1 implies that 1/cj = hj — aj. 


Note that the types of densities may not respect the order of 
sensors in [0, L], e.g. the 1st and 3rd distances can be from the 
2nd group of densities equal to /2(/), while the 2nd distance 
can be from the 1st group. In this case we say that index 1 
belongs to group 2, symbolically (1) = 2. Here the brackets 
(•) denote the operator transforming an index ? = 1, . . . , n of 
a distance into its group number (i) varying from 1 to k. 

For a heterogeneous network, the function Vn{r,l) from 
section III will be a sum over arrays of signs Q = 
{qi,. . . ,qn) depending on prescribed densities {/i, . . . , /fc}. 
Let [aj{r),bj{r)] be the intersection of [0,r] with [aj,bj], 
where fj ^ 0. Set qf = (1 ± qi)/2, e.g. 1+ = 1, 1" = 0. 

The Heterogeneous Corollary. In the above notations and 
under the conditions of the Connectivity Theorem assume 
that distances between successive sensors have probability 
densities fj{l) = Cj on [aj,bj], j — Then the 

Vn{R, L) 



probability of cormectivity is P„ = 



Vn{L,L) 



, where 



{Q)<l 

<9=(9l,-,«n) 



{1 - {Q)Y and 



i=l 



i=l 



The indices in the brackets {i) from the last formula above 
take values 1, . . . , A: for each i = 1, . . . ,n, i.e. [a(i) (r), (r)] 
is the segment where the («)-th density /(j) is defined after 
restricting it to the transmission range [0,r]. In particular, if 
each i-th distance has its own density then {i) = i and the 
indices i, j = 1, . . . ,n are equal to each other. 

First we show that the Constant Corollary is a very partial 
case of the Heterogeneous Corollary with only one constant 
density /i = l/{b — a) on [a,b], i.e. fc = 1. To compute 
Vn{L,L) we note that [ai{L),bi{L)] = [a,b]. Let k be the 
number of pluses in an array Q. Then dq = {—l)'^/{b — 
a) and (Q) = a{n — fc) + bk. So the sum over Q can be 
rewritten as a sum over < fc < n. For any fixed fc there 
are (^) different arrays Q containing exactly k pluses. By 
the Heterogeneous Corollary the common term in the sum 
Vn{L,L) over fc is {-!)'' {1){L - a{n - k) - bk)". The only 
difference in computing Vn{R,L) is that 6i(i?) = R, which 
leads to the formula from the Constant Corollary. 



The complexity to compute the function Vn{r,l) from the 
Heterogeneous Corollary is 0(2"), because Vn{r,l) is a sum 
over 2" arrays of signs and (Q) is a weighted sum of 
endpoints ai{r),bi{r). In the general case, the expression (Q) 
can take 2" different values. If there are only fc different 
endpoints then the algorithm has the polynomial complexity 
0{n''), see the 3-step Density Corollary in Appendix D. If all 
the segments [aj,bj] are subsets of [0,R] then any network 
will be connected and the formula above gives 1, because 
the numerator of P„ coincides with the denominator when 
aj{R) = aj{L) and bj{R) = bj{L), j = l,...,k. 

Proof of the Heterogeneous Corollary extends the proof of the 
Constant Corollary. We consider the truncated densities 

= C(,)(^i(? - a(,)(r)) - u{l - 6(,)(r))), 

i = 1, . . . , n, where (j) denotes the group containing the ith 
distance. By Lemma 2(c) for m = one has 

LT{/M}(a) = c(,)^ ^ . 

Substitute each Laplace transform LT{/|'^'}(s) into the func- 
tion g{s) from Lemma 3 and expand the product g{s), which 
gives the following sum of 2" terms: 



9{s) 



-A[ ^ = 2^ '^Q - 



„ri+l 



i=l 



Q 



The sum is taken over arrays Q = {qi,...,qn) of signs. 
The sign = — 1 means that the term with a(j)(r) is 
taken from the i-th factor, the sign qi = +1 encodes the 
second term with &(i)(r). The total power of the exponent 
in the resulting term corresponding to Q is —{Q)s, where 

n 

(Q) = J2{(^(i){^)'li + b(i){r)q^). So each minus contributes 

i=l 

— a(j)(r)s to the total power, while each plus contributes 
—6(j) (r)s. Each plus contributes factor (—1) to the coefficient 

dq, i.e. = H as required. 

i=l 

Compute the inverse Laplace transform by Lenuna 2(d): 

vn{r, I) = LT-\g{s)} = ^ ^(^ - (g))"^^ - (0)), 

where the unit step functions u{l — {Q)) can be replaced by 
the upper bound I < {Q) as in the final formula. □ 

The algorithm for computing the function w„(r, I) from the 
Heterogeneous Corollary has the following steps: 

• initialise 2 arrays a(j)(r), (r), i = 1, . . . ,n; 

• make a computational loop over 2" arrays Q of signs; 

• for each Q compute (Q) and check the upper bound I < {Q), 
find dQ{l — (Q))" and add it to the current value of Vn{r, I). 

The algorithm for computing Vn{L, L) is similar, replace R 
by L. If we are interested only in P„, we may forget about n! 
which is canceled after dividing Vn{R, L) by Vn{L, L). 
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Appendix D 
PiECEWisE Constant Densities 

In this appendix we show how to compute the probability 
of connectivity building any piecewise constant density from 
elementary blocks in the Heterogeneous Corollary. The build- 
ing engine is the Average Density Corollary below dealing 

k 

with the average /(s) — ^ fj(s)/k of constant densities 

fj{s) = Cj on [aj,bj] and fj{s) = otherwise, j = 1, . . . ,k. 
The factor 1/k guarantees the condition f{s)ds = 1, which 
follows from fj{s)ds = 1, j = 1, . . . ,k. 

For any ordered partition n = ni + ■ ■ ■ + Uk into k non- 
negative integers, denote by {rii, . . . ,nk) the collection of 
densities, where the first m densities equal /i, the next n2 
densities equal /2 etc. For example, given 2 constant densities 
/i, /2, number n = 3 can be split into 2 non-negative integers 
in one of the 4 ways: 3 = 0-1-3=1 + 2 = 2 + 1= 3-^0. 
Then (1,2) denotes the collection (/i,/2,/2), i-C- the 1st 
distance in such a network has the density /i, while the 
remaining 2 distances have the density /2. For each partition 
{rii, rifc) or, equivalently, a collection of constant densities, 
let ^(^"1 ' ■ (^-f, I) be the function defined by the formula from 
the Heterogeneous Corollary in Appendix C. 

The Average Density Corollary. In the above notations and 
under the conditions of the Connectivity Theorem if distances 
between successive sensors have the probabihty density f{l) = 

k 

J2 fj{l)/k on [0, L], then the probabihty of connectivity is 

E^i""-'"''^(i^,i^)M!...nfc! 
Pn = 7 ^ . Both sums are taken 

over all collections of densities (ni, . . . , n^) corresponding to 
ordered partitions n = ni + • • • + n^. 

The products ni! . . . n^! can not be canceled in the formula 
above, because the numerator and denominator of Pn are sums 
of many terms involving different products ni\...nk\ over 
all ordered partitions n = ni + • • • + rife. The complexity 
to compute P„ is 0(n2"), because each function 
is computed by the algorithm describe after the Heteroge- 
neous Corollary using 0(2") operations. In partial cases 
the computational complexity can be reduced to polynomial, 
see comments after the 3-step Density Corollary below. The 
algorithm computing the probabihty from the Average Density 
Corollary applies the algorithm from the Heterogeneous Corol- 
lary to each function L) and vt^'-'""'\L, L) 
substituting the results into the final formula above. 

Proof of the Average Density Corollary. 

We may forget about the factor 1/fc as usual. Set gj{s) = 

LT{/|'^'}, j = l,...,k. Lemma 3 impUes that Vn{r,l) = 

k 

LT~^{g'(s)}, where ^(s) = {Y^ 5j(s))"/s. Expand the brack- 



ets: g{s) = J2 



n\\ . . .rife! 



, where the sum is over 



all partitions n = n\^ h rife into k non-negative integers. 

By Lemma 3 each term g"^ . . . /s is the inverse Laplace 
transform of the function w^"^' 0' where the first n\ 



distributions are /i, the next distributions are etc. It 
remains to cancel n\ in the final expression. □ 

By taking sums of constant densities Cj on [a^, hf\, one can 
get any piecewise constant function on [0, L]. Any reasonable 
function can be approximated by piecewise constant ones. 
Hence the Heterogeneous Corollary and Average Density 
Corollary are building blocks for computing the probabihty 
of connectivity for any real-life deployment of sensors. 

We demonstrate this universal approach for the sum of 2 
constant densities over 2 different segments. So the density 
in question is a 3-step function depending on the radius R 
and one more parameter C, its graph is shown in Fig. 9. Let 
f = {f^+ /2)/2 be the density on [0, L] such that 

C ifle[0,R], 
otherwise; 



A (0/2 



f2il)/2- 



1/R-C if I e [R/2,3R/2], 
otherwise. 



where 0<C<l/i?isa constant and 3i?/2 < L, see Fig. 9. 

L 

The constants C and 1/i? are chosen so that / f{l)dl = 1. 





/(O 








1/R 








1/R-C 




c 








R/2 


R 


3R/2 


I 



Fig. 9. The piecewise constant distribution depending on R, C 

From Fig. 9 for a network of a sink node at and 1 sensor 
at yi the probability of connectivity is P(0 < yi < R) = 
(CR + l)/2, the area of the first two rectangles below the 
graph of f{l). For example, if C = 0.9/R then Pi = 0.95 
as shown in Fig. 10, so it is very likely that 1 sensor will 
be close enough to the sink, although such a network can not 
cover the whole segment [0, L]. The 3-step Density Corollary 
below gives an example how to compute the probability of 
connectivity explicitly for a piecewise constant density using 
the Heterogeneous Corollary and Average Density Corollary. 

The 3-step Density Corollary. Under the conditions of the 
Connectivity Theorem and for the piecewise constant density 
f{l) above, the probability of connectivity is P„ = 

n m n-m (_l)fei+fc2 _ (^2ki + ^2 + 71 - TO)i?/2)" 

m=o fei=ofe2=o dmki\{m - fci)!fc2!(n - m - ^2)! 



rii n—m 



E E E 

m=0/ci=0 fe2=0 



(■_l)fei+fe2 _ ^2ki + 2fc2 + n - m)R/2Y ' 
dmk\\{m — fci)!fc2!(n — m — A;2)! 
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where dm = C~™{l/R — C)™^", the sums are over all 
possible values of m,ki,k2 such that the expressions in the 
brackets taken to the power n are positive. 

The complexity to compute the probability P„ above is 
O(rV^), because the sums in the numerator and denominator 
are over 3 non-negative integers not greater than n and each 
term requires 0(1) operations. If C ~ l/R, i.e all distances 
are in [0, i?], then set dm ~ for m < n. Hence m — n, 
k2 = and the sums over 3 parameters m,ki,k2 reduce to 
the same single sum over fci = 0, . . . , n in the numerator and 
denominator, which gives P„ = 1 as expected for C = l/R. 

If C = 0, i.e. each distance is uniformly distributed on 
[R/2, 3i?/2], then set = for to > 0. Therefore, to = 0, 
ki — and the result containing only sums over k2 = 0, . . . ,n 
coincides with the probability P^l from the Constant Corollary 
with [a,b] — [R/2,3R/2] after canceling do and multiplying 
the numerator and denominator by ^2! to get (^^) 

n 

S (-1)'=^ {L - {k2 + n)R/2)^/k2\{n - ^2)! 
J2 - (2fc2 + n)R/2)^/k2l{n - fca)! 

k2=0 

In the 3-step Density Corollary for n = 1 both sums contain 
only 4 non-zero terms corresponding to the parameters 

(to, fci, fca) = (0,0,0); (0,0,1); (1,0,0); (0,1,0). 

Then all the factorials in the formula are 1 and we get Pi ~ 
(^ - C)(£ - f ) - (^ - C){L -R) + CL~ CjL - R) 
(i-C)(L-f)-(^-C)(L-f ) + CL-C(L-i?) 

— {CR + l)/2 as we have checked using Fig. 9 directly. 

Proof of the 3-step Density Corollary. 

In the notations of the Heterogeneous Corollary we have only 
k = 2 densities. Let the first to distances between successive 
sensors have the probability density /i, while the last n — m 
distances have the density /2. An array Q of n signs similarly 
splits into two parts consisting of to signs and n — m signs. 
Let ki and /c2 be the number of pluses in each part. 

To compute ii,^™'""™-' [L, L) from the Heterogeneous Corol- 
lary for the partition n = m + {n — m) we note that 

[ai(L),6i(L)] = [0,R], [a2(L), &2(i)] = [i?/2, 3i?/2], 

dg = /dm, (Q) = (fci + k2 + {n- m)/2)R. 

So the sum over arrays Q can be rewritten as a sum over ki,k2. 
For fixed values of these parameters, there are 
different arrays of signs. After canceling the factorials to! 
and (n — m)! in the Average Density Corollary, the sum 

J2 " ™ {L, L) is the required denominator of P„. 

m=0 

The only difference in computing Vn{R, L) is that 62 (i?) = 
R, not 3i?/2. This replaces 2k2 by k2 in the numerator. □ 

Given the piecewise constant density f{l) with the inter- 
mediate parameter C — 0.9/ R, Table 6 shows the minimal 
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Fig. 10. The probability of connectivity for f{l) with C = 0.9/i? 

number of sensors having different radii such that the network 
in [0, L] is connected with probability 0.95, where L = 1km. 
Fig. 10 shows the probability of connectivity P„ > 0.2. 

Table 6. The probability of connectivity for the piecewise 
constant density with C = 0.9/i? and different radii. 



Transmission Radius, m. 


250 


200 


150 


100 


50 


Min Number of Sensors 


12 


17 


25 


44 


105 
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